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^ ■ Abstract 

> : We report about an ana ly fe study involving the interm^te wave packet formalism for q uau- 

tifying the physically relevant information which appear in the neutrino two-flavor conversion 
■ formula and help us to obtain more precise limits and ranges for neutrino flavor oscillation. By 

following the sequence of analytic approximations where we assume a strictly peaked momentum 
distribution and consider the second-order corrections in a power series expansion of the energy, 

o 

we point out a residual time-dependent phase which, coupled with the spreading /slippage effects, 



(N 
O 



can subtly modify the neutrino oscillation parameters and limits. Such second-order effects are 
usually ignored in the relativistic wave packet treatment, but they present an evident dependence 
Q_i| on the propagation regime so that some small modifications to the oscillation pattern, even in the 



Oh 



ultra-relativistic limit, can be quantified. These modifications are implemented in the confront 
with the neutrino oscillation parameter range (mass-squared difference Am 2 and the mixing-angle 



9) where we assume the same wave packet parameters previously noticed in the literature in a 
kind of toy model for some reactor experiments. Generically speaking, our analysis parallels the 
recent experimental purposes which concern with higher precision parameter measurements. To 
summarize, we show that the effectiveness of a more accurate determination of Am 2 and 9 depends 
on the wave packet width a and on the averaged propagating energy flux E which still correspond 
to open variables for some classes of experiments. 
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I. INTRODUCTION 



In the last years, the quantum mechanics of neutrino oscillations [l, 2, 3] has experienced 
much progress on the theoretical front [41, not only in phenomenological pursuit of a more 
reflned flavor conversion formula fl fl Q, but also in efforts to g ive the theory a fo rI „a. 
structure within quantum field formalism ay , [hJ. From the point of view of a first 
quantized theory and in the context of vacuum oscillations, as a first analysis, the probability 
that neutrinos originally created as a u a flavor-eigenstate with averaged energy E oscillates 
into a up flavor-eigenstate over a distance L is given by 

P(y a — > up; L) = sin 2 (20) sin 

where we have assumed that the main aspects of the oscillation phenomena can be under- 
stood by studying the simple-minded two-flavor problem constructed in terms of the u x (x) 
and u 2 (x) plane- wave mass-eigenstates with the two-flavor mixing angle represented by 6 and 
the mass-squared difference given by Am 2 .. As appointed by the D. Groom PDG review 111 ]. 



Ami L 
4E 



, (c,H=l) 



(1) 



although this equation is frequently quoted and used in Monte Carlo calculations, the wave 
function is badly behaved for reasons larger than about one, where it oscillates more and 
more rapidly in the interval between sin (26) = and sin (26) = (P) as the argument 



4£ 

increases. Moreover, it is difficult to relate this function to the exclusion curves described 
in the literature 



111- 



In fact, the intermediate wave packet (WP) approach |12J eliminates the most controver- 
sial points rising up with the standard plane- wave (PW) formalism [l3,[l4| since wave packets 
describing propagating mass-eigenstates guarantee the existence of a coherence length 12], 



avoid the ambiguous approximations in the PW derivation of the phase difference 



151 ] and, 



under particular conditions of minimal slippage between the mass-eigenstate wave packets, 
recovers the oscillation probability given by the standard PW treatment 39]. Enclosed by a 
restrict motivation, we avoid a more extensive discussion about several controversial aspects 
concerning with the intermediate WP formalism (see for instance 4j, ll7| which work on the 
external WP framework) and, quite generally, we follow an analytical approach where the 



mass-eigenstate time evolution does not concern wit 



i the WP limitations. As preliminary 



investigation we consider gaussian wave packets , as we know, the unique which 

enable us to analytically quantify the first and the second-order corrections to the flavor 
conversion formulas. By means of the second-order terms we can compute the effects of 



an extra time-dependent phase which is added to the standard oscillation term 



Am 2 t 
2E 



and modifies the oscillating character of the propagating particles. Generically speaking, we 
avoid more sophisticated methods as field theoretical prescriptions {4, 18] in detriment to 
a clearer treatment which commonly simplifies the understanding of physical aspects going 
with the oscillation phenomena. 

Reporting to such an analytical construction [19], we look for improving the procedure 
for obtaining the two-flavor oscillation parameter exclusion region boundary for a generic 
class of oscillation experiments. In fact, precise determination of the oscillation parame- 
ters and search for non-standard physics such as a small admixture of a sterile component 
in the solar neutrino- flux are still of interest. In addition, from the experimental point of 
view, to determine Am 2 more precisely, further KamLAND exposure to the reactor neu- 



trinos are becoming most powerful 



20 



211 ] at the same time that more precise neutral 



current measurements by SNO can contribute in reducing the uncertainty of the mixing 



angle 



21 



22] . By knowing the emergence of neutrino physics is fueled by such a recent 



growth in quality and quantity of experimental data, our main purpose is discussing how 
much the determination of mixing parameters and mass- differences can be improved in the 
context of the wave packet phenomenological analysis and, eventually, suggest a perspective 
of improvement on bounds on 9 13 before experiments designed specifically for this parameter 
start. The first step of our study, which is presented in section II, concerns with the an- 
alytical derivation of a flavor conversion formula where a gaussian momentum distribution 
and a power series expansion of the energy up to the second-order terms are introduced 
for obtaining analytically integrable probabilities. Adopting a strictly peaked momentum 
distribution to construct each mass-eigenstate wave packet allows us to analytically quan- 
tify these rising-up second-order corrections. Firstly, we reproduce the localization effects 
described in terms of the well-known wave packet spreading and of the slippage between the 
wave packets, which leads to the decoherence between the propagating mass-eigenstates. In 
particular, it is well-established that the decoherence effect has a correspondence with an 
oscillation damping (exponential) factor 12J, |23j, thus we just report about such effects to 
show that both of them can be quantified, up to second-order corrections, in non-relativistic 
(NR) and ultra-relativistic (UR) propagation regimes. In parallel, we also recover an addi- 
tional time-dependent phase which changes the standard oscillating character of the flavor 
conversion formula. However, the main contribution of our study, which is presented in 
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section III, concerns with the understanding of the physical aspects carried by neutrino 
two-flavor oscillation parameters and limits. In particular, we suggest the existence of more 
precise corrections in the procedure 11( where the wave packet characteristics are used for 
describing the mass-squared difference Am 2 and the mixing-angle 9 which are constrained 
by the experimental input parameters. In section IV, by the time we assume the same test- 



ing wave packet parameters previously introduced in the literature 11( (for some reactor 
experiments), we can establish a phenomenological comparison which allows us to confront 
the energy second-order corrections introduced in the previous sections with the predicted 
PW and WP with first-order correction results. We also observe the possibilities for ex- 
tending the analysis to other production- type neutrinos (solar and supernova) . Finally, we 
draw our conclusions in section V by emphasizing that the effectiveness of a more accurate 
determination of Am 2 and 9 can eventually depend on the wave packet width a and on the 
averaged propagating energy flux E. 



II. WAVE PACKETS WITH SECOND-ORDER CORRECTIONS 



The time evolution of flavor wave packets can be described by the state vector 

&(z,t) = cf) 1 (z,t) cos 9 ^ + <p 2 (z,t) sin 6 u 2 

= [0i (z, t) cos 2 9 + 4> 2 {z, t) sin 2 9] u a + [0 2 (z, t) — (p 1 (z,t)] cos 9 sin 9 isp 

= (j) a (z,t;9)v cl + 4) !3 (z,t;9)vf3, (2) 

with flavor and mass-eigenstate indices as previously defined. The probability of finding a 
flavor state Vp at the instant t is equal to the integrated squared modulus of the Vp coefficient 



P{v a -> Up] t) 



dz |^(z, t;9)\ 2 = { i _ lNT(t) } ? 



where iNT(t) represents the mass-eigenstate interference term given by 



iNT(t) = Re 



dz 4>l(z,t) (p 2 (z,t) 



Let us consider mass-eigenstate wave packets given by 



<M^°) = ( —A ex P 
1 7rcr / 



exp [ipi z] 



(3) 



(4) 



(5) 



at time t = 0, where i — 1, 2. The wave functions which describe their time evolution are 



<t>i(z,t) 



dp z 
2ir 



- iE p ] J + 



ip z z 



(6) 



where Ep] = (p 2 z + m 2 ) 2 and if(j) z -p i ) = (27ra 2 ) 4 exp 



(Pz-Pi) 2 a 2 



. In order to obtain the 



oscillation probability, we can calculate the interference term Int(£) by solving the following 
integral 



/ ~2jT ^ Pz ' Pl ^ Pz ~ P2 ^ eXP H AE Pz *] 



exp 



-(aAp) 2 



— ^ ip 2 {p z -Po) exp [-i AE Pz t], 



(7) 



where we have changed the ^-integration into a p z -integration and introduced the quantities 
Ap = p 1 -p 2 - l p = \{Pi + P2) and AE Pz = Ep 1 } -EpJ . The oscillation term is bounded by 
the exponential function of aAp at any instant of time. Under this condition we could 
never observe a pure flavor-eigenstate. Besides, oscillations are considerably suppressed if 
aAp > 1. A necessary condition to observe oscillations is that aAp <C 1. This constraint 
can also be expressed by 5p 3> Ap where 5p is the momentum uncertainty of the particle. 
The overlap between the momentum distributions is indeed relevant only for bp ^> Ap. 
Strictly speaking, we are assuming that the oscillation length (-7T , 3 ) is sufficiently larger 
than the wave packet width which simply says that the wave packet must not extend as 
wide as the oscillation length, otherwise the oscillations are washed out [12I . 
Turning back to the Eq. ([7]), without loss of generality, we can assume 
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lNT(t) = Re {j ^ f 2 {p z -p ) exp [-» AE Pz t] j 



(8) 



In the literature, this equation is often obtained by assuming two mass-eigenstate wave 
packets described by the same momentum distribution centered around the average momen- 
tum p = p . This simplifying hypothesis also guarantees instantaneous creation of a pure 
flavor eigenstate u a at t = [15J. In fact, for (fi^z, 0) = (p 2 (z, 0) we get from Eq. ([2]) 



( 2 V 


' z 2 ' 


— 2" ex p 




\7ra 2 J 


~ a 2 



exp [ip z] 

and <fii3(z,0,8) = 0. In order to obtain an expression for (p^Zjt) by analytically solving the 

(0 



integral in Eq. ([6]) we firstly rewrite the energy EpJ as 



-(0 



Ei 



1 + 



Pt-Po 
E 2 



(9) 



Pz-PQ 

Ei ■ 



where Ef = (m 2 + pi) 2 , v< = fr and a, 

We are attentive to the fact that the assumption of wave packets with gaussian shape can 
introduce some (not quantifiable) limitations to the interpretation of the following results. 



However, in this kind of analysis and, for instance, in non-relativistic quantum mechanics, 
free-propagating gaussian wave packets are frequently assumed because the calculations can 
be carried out exactly for these particular functions and, consequently, the main physical 
aspects can be easily interpreted from the final analytically obtained expressions. The reason 
lies in the fact that the frequency components of the mass-eigenstate wave packets, Ep} = 
p 2 z /2m i , modify the momentum distribution into "generalized" gaussian, easily integrated by 
well-known methods of analysis. The term p 2 z in Ep} is then responsible for the variation in 
time of the width of the mass-eigenstate wave packets, the so-called spreading phenomenon. 
In relativistic quantum mechanics the frequency components of the mass-eigenstate wave 
packets, Ep] = ^p 2 + m 2 , do not permit an immediate analytic integration. This difficulty, 
however, may be remedied by assuming a sharply peaked momentum distribution, i. e. 
(aE^ 1 ~ <Ji <C 1. Meanwhile, the integral in Eq. ([6]) can be analytically solved only if 
we consider terms up to order af in the series expansion. In this case, we can conveniently 
truncate the power series 



E, 



07 



1 + OiV, + -f (1 



m 



(10) 



and get an analytic expression for the oscillation probability. The zero-order term in the 
previous expansion, E t , gives the standard plane- wave oscillation phase. The first-order 
term, p^a^ is responsible for the slippage between the mass-eigenstate wave packets due to 
their different group velocities. It represents a linear correction to the standard oscillation 

2 

phase 15]. Finally, the second-order term, , which is a (quadratic) secondary correction 
will give the well-known spreading effects in the time propagation of the wave packet and 
will be also responsible for an additional phase to be computed in the final calculation. In 
case of gaussian momentum distributions, all these terms can be analytically quantified. By 
substituting (ITU]) in Eq. (jHJ) and changing the ^-integration into a ^-integration, we obtain 
the explicit form of the mass-eigenstate wave packet time evolution, 



<t> t (z,t) 



where z) 



arctan 



7TO?(t) 

2 m 2 t 



exp [-i (9i(t, z) + Et t - p z)} exp 
and ai(t) 



a 2 (t) 



(11) 



2m\t (z-Vj t) 2 



a 1 + 



t 2 J 2 . The time- 



dependent quantities a^t) and O^t, z) contain all the physically significant information which 



arise from the second-order term in the power series expansion ({TO]) . By solving the integral 
(JHJ) with the approximation (Q and performing some mathematical manipulations, we obtain 



iNT(t) = BND(i) X OSC(t), 



where we have factored the time-vanishing bound of the interference term given by 
and the time-oscillating character of the flavor conversion formula given by 



(12) 



(13) 



Osc(t) = Re{exp[-iAEt-ie(t)}} 
= cos[AEt + 0(t)], 



where 



and 



S p( t ) = -A ( ^- 



Avt 

a 2 p 



©(*) 



\ arctan [Sp(<)] 



a 2 p 2 



SP 3 (t) 



V [l+S P 2 (i) 



(14) 



(15) 



(16) 



with p = 1 - [3 + (^) 2 ] ffc, AE = £?W - E& and E = y/E~E~ 2 . The time-dependent 
quantities Sp(£) and Q(t) carry the second-order corrections and, consequently, the spreading 
effect to the oscillation probability formula. If AE -C E, the parameter p is limited by the 

2 

interval [1,-2] and it assumes the zero value when ^ ~ |- Therefore, by considering 
increasing values of p , from non-relativistic (NR) to ultra- relativistic (UR) propagation 
regimes, and fixing J^ 2 , the time derivatives of Sp(£) and G(t) have their signals inverted 

2 

when -|| reaches the value |. The slippage between the mass-eigenstate wave packets is 
quantified by the vanishing behavior of Bnd(£). In order to compare Bnd(£) with the 
correspondent function without the second-order corrections (without spreading) , 

(Avt) 2 " 



BND ws (t) = exp 



2a 2 



we substitute Sp(t) given by the expression (fill) in Eq. ([IB"]) and we obtain the ratio 



BND(t) 
BND ws (t) 



exp 



2 a 6 E s 



2 AEt 



(17) 



(18) 



The NR limit is obtained by setting p 2 = 1 and p = in Eq. ( jT7|) . In the same way, the 
UR limit is obtained by setting p 2 = 4 and p = E. In fact, the minimal influence due to 



second-order corrections occurs when fa 



|(P«0). 
7 



The oscillating function Osc(i) of the interference term Int(£) differs from the standard 
oscillating term, cos [AEt], by the presence of the additional phase Q(t) which is essentially 
a second-order correction. The modifications introduced by the additional phase G(t) are 
discussed in Fig. [1] where we have compared the time-behavior of Osc(i) with cos [AEt] 
for different propagation regimes. The effective bound value assumed by 0(t) is determined 
by the vanishing behavior of Bnd(£). To illustrate this flavor oscillation behavior, we plot 
Doth the curves representing Bnd(£) and Q(t) in Fig. [21 (both figures are obtained from 
3]). We note the phase slowly changing in the NR regime. The modulus of the phase |0(£)| 
rapidly reaches its upper limit when || > ~ and, after a certain time, it continues to evolve 
approximately linearly in time. Essentially, the oscillations vanishes rapidly. By superposing 
the effects of Bnd(£) in Fig. [2jand the oscillating character Osc(t) expressed in Fig. [TJ we 
immediately obtain the flavor oscillation probability which is explicitly given by 

P{y a - v# t) « |l - [1 + Sp 2 (t)H exp \- 2a 2 ^f p2 ( J cos [AEt + &(t)}\ (19) 

and illustrated in Fig. [31 Obviously, the larger is the value of a E, the smaller are the wave 
packet effects. 

III. UNDERSTANDING TWO-FLAVOR OSCILLATION PARAMETERS AND 
LIMITS 

As earlier extensively discussed in the literature 0, [5, I2I, 16], there are various restrictive 
conditions under which the two-neutrino mixing approximation is valid and the standard PW 
as well as several classes of WP treatment can be utilized for describing the flavor conversion 
phenomena. Respecting such limitations, in the previous section we have quantified the 
second-order corrections introduced by an intermediate WP analysis. Now, we are in position 
to verify how such modifications can affect the reading of the neutrino oscillation parameter 
ranges excluded by the confront with experimental data. Therefore, our effective purpose 
is contrasting the standard PW treatment with the WP analytic study presented in the 
previous section by re-obtaining the two-flavor neutrino oscillation parameters and limits in 
a very particular phenomeno logical context. 

From a practical point of view, we have to establish the input experimental parameters 
as being the detection distance L (from source), the neutrino energy distribution E and the 
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appearance (disappearance) probability (P). In addition, to make clear the initial propo- 
sition, it is instructive to redefine some parameters which shall carry the main physical 
information in the oscillation formula, i. e. 

1 L 



b 



A. 

aE 



and v 



Po 
E ' 



(20) 



2v E 1 + E 2 

with E previously defined. For real experiments, E and L can have some spread due to 
various effects, but in a subset of these experiments, there is a well-defined value of b about 
which the events distribute ll|]. Following the same approach we have adopted while we 



were analyzing the parameter p in Eq. (1151) . if AE <C E, which is perfectly acceptable from 
the experimental point of view, we can write E = y/E x E 2 ~ \(E\ + E 2 ) so that an effective 
PW flavor conversion formula can be obtained from Eq. (TTJ) as 

sin 2 (29) 



P(u a -> vp;b ) 



{1 - cos [2b Am 2 ]} 



(21) 



where b carries the dependence on the detection distance L and on the propagation regime 
(v). Performing some analogous substitutions, the interference terms of Eq. (fl2l) which 
explicitly appears in the WP flavor conversion formula of Eq. (1191) can be read in terms of 
the above rewritten parameters 

28 2 b y2 ( Am ") 



Bnd(6 ) = [l + Sp 2 (& )Hexp 



1 + Sp 2 (b 



and 



Osc(6 ) = cos[2 6 Am 2 + e(6 )], 



where 



and 



SI 



e(6 ) 



Sp(6 ) = -^Am 2 

"0 



1 ro /, m blv 2 Sp 3 (6 ) 

- arctan [S P (6 )]- 2 ^ p2[1 + Sp2(M] 



(22) 



(23) 



(24) 



(25) 



with p m 1 — 3v 2 . Attempting to the rate a = Sf,/b = (aE)' 1 which carries the relevant 
information concerning with the wave packet width and the averaged energy flux, if it is was 
sufficiently small (a 1) so that we could ignore the second-order corrections of Eq. (J9j), 
the probability only with the leading terms could be read as 

sin 2 (29) 



{1 - cos [2b Am 2 } exp [-2(5 b v Am 2 ) 2 ]} 



(26) 



which brings up the idea of a coherence length L co h ~ [4], |23| and, in the particular case 



of an UR propagation (v = 1), is taken as a reference in the confront with experimental 
data By the way, despite the relevant dependence on the propagation regime (v), once 
we are interested in some realistic physical situations, the following analysis will be limited 
to the UR propagation regime corresponding to the effective neutrino energy of the current 
flavor oscillation experiments. 

Strictly speaking, most results in the neutrino mixing listings are presented as Am 2 limits 
(or ranges) for sin 2 (29) = 1 , and sin 2 (26 l ) limits (or ranges) for large Am 2 . Together, they 
summarize the most of the information contained in the usual Am 2 x sin 2 (29) plots which 
provide the parameter exclusion region boundary in the experiments' papers. Thus, we 
can compare the PW and WP resolutions by enumerating some relevant aspects which can 
be observed from the curve Am 2 x sin 2 (29) in the Fig. H] where the Palo Verde 25| and 
the KamLAND [26j experiments/exclusion curves are represented. For the moment, we are 
interested in the analytical characteristics of such exclusion curves, which, however, we turn 
to analyze under the phenomenological point of view in section IV. The main aspects to be 
quantified are: 

1) In both PW and WP cases, for large Am 2 the fast oscillations are completely washed 
out respectively by the plane-wave resolution or by the smearing out behavior due to the 
mass-eigenstate wave packet decoherence. Consequently sin 2 (29) = 2(P) in this limit. 

2) For PW calculations the maximum excursion of the curve to the left occurs at 
Am 2 = 7r/26 when sin 2 (2$) = (P). When the WP Eq. ( |26l) is used such a maximal point 
occurs at the solution of the transcendental equation — 2cr 2 Am 2 b = tan [2 Am 2 b ] which 
can be approximately given by Am 2 « n/2b (1 + 0(a 2 )), but we know that the second- 
order terms (a 2 ) are not being considered. If we had taken into account the second-order 
corrections in the WP analysis, the maximal value would have been accurately given by 
Am 2 « Tc/2b (1 + 2cr 2 + 0(a 4 )), consequently, a little smaller value than the PW solution. 

3) By qualitatively assuming the well-established phenomenological constraints which set 
Am 2 <C leV 2 when sin 2 (29) ~ 1, we can reconstruct the nearly straight-line segment at 
the bottom of the curve by expanding the probability expressions up to order 0((Am 2 ) 2 ) so 
that we can obtain the generic solution 



Am 2 « V( P ) ^ (27) 

b sm(26)^F App (a) 
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where we have F PW (cr) = 1 for the PW limit, F WP1 (cr) = 1 + o 2 for the WP treatment with 
first-order corrections and F WP2 (a) = 1 + 2<t 2 + |cr 4 for the WP treatment with second-order 
corrections. Eventually, if one had abandoned the analytic calculations and had taken into 
account higher order terms in the power series expansion of Eq. ( fTOl) . there would have 
been some minor corrections to the a 4 term in F WP2 - ^Ve discard such minor corrections 
by assuming F WP2 (cr) ~ 1 + 2cr 2 which, in fact, is the correct approximation when we are 
considering the energy expansion up to second-order terms and sufficient for comparing the 
approximations in the Fig. [51 We emphasize that we must consider the second-order term 
in the series expansion in the Eq. ( ITUl) since the modifications emerge with a 2 . As a testing 
case, we consider the toy model presented in llj where the values of a = 0.23 and a = 0.3 
are used for plotting the curve Am 2 x sin 2 (29) . From an immediate analysis of the Eq. ( |27j) 
with WP2 approximation we can conclude that more accurate values of Am 2 are constrict 
to be diminished by approximately 8% of the value computed with the PW approximation. 

Quite generally, the complete analysis of the amortisation coupled to the oscillating char- 
acter depends upon the experimental features such as the size of the source, which allows us 
estimating the wave packet width (a), the neutrino energy distribution (E) and the detector 
resolution (L ). Anyway, the main point to be considered is that the characterization of the 
wave packed (a) implicitly described by a, which is accompanied by the precise determina- 
tion of the neutrino energy distribution (E), plays a fundamental role when the neutrino 
oscillation parameter and limit accuracy is the subject of the phenomenological analysis. To 
clarify this point, in addition to the information interpreted from Fig. [5], in the next section 
we turn back to the phenomenological discussion on the Fig. H] which illustrates one simple 
case of parameter determination for reactor neutrinos by considering the results from some 



disappearance experiments 



25 



26]. 



IV. PHENOMENOLOGICAL CONSTRAINTS TO REACTOR EXPERIMENTS 
EXTENSION TO SOLAR AND SUPERNOVA NEUTRINOS 



The first hints that neutrino oscillations actually occur were serendipitously obtained 
through earl y st udies of solar neutrinos 271] and neutrinos produced in the atmosphere by 



cosmic rays 30j, |3l|. More recently, nuclear reactors and particle accelerators have consti- 



tuted another source of neutrinos utilized for accurate measurements of flavor oscillation 
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parameters and limits. Reactor neutrino experiments correspond essentially to an electron- 
antineutrino v e disappearance experiment where, generically speaking, one looks for the 
attenuation of the initial neutrino flavor-eigenstate u a beam in transit to a detector, where 
the u a is measured. In contrast to the detection of even a few wrong-flavor (i/g) neutrinos es- 
tablishing mixing in an appearance experiment, the disappearance of a few right-flavor (u a ) 
neutrinos goes unobserved because of statistical fluctuations For this reason, disap- 

pearance experiments usually cannot establish small-probabilities (sin 2 (29) <C 1). Besides, 
they can fall into several situations llj into which we do not intend to go deep. 

By following the purpose of a comparative phenomenological study, we take into account 
the experimental data from the Palo Verde 25j] and KamLAND [26] experiments by means 



of which exc 
periment 



usion plots in the plane Am 2 x sin 2 (29) can be elaborated. The Palo Verde ex- 



25] consists in a disappearance search for v e oscillations at 0.75 — 0.89 km distance 



from the Palo Verde reactors. As consequence of the experimental analysis we assume the 
input parameter (P) < 0.084 and the averaged energy E set in the interval 3.5 — 4.2 MeV. 
The KamLAND collaboration [26] observes reactor v e disappearance at ~ 180 km baseline 
to various Japanese nuclear power reactors. This is the lower limit on the mass difference 
spread unlike all other disappearance experiments [26] and the observation is consistent with 
neutrino oscillations, with mass- difference and mixing angle parameters in the Large Mixing 
Angle Solution region of the solar neutrino problem. In this case, we assume the input pa- 
rameter (P) > 0.2 with the reactor u e energy spectrum smaller than 8 MeV and an analysis 
ihreshold of 2.6 Me V where the experiment sensitive Am 2 range is set down to 10~ 5 eV' 2 
261 ] . The important point we attempt in the Fig. H] is the comparative modifications due to 
second-order corrections, i. e. we are not setting extremely accurate input (experimental) 
parameters but we are setting a more accurate procedure for obtaining the output parameter 
(Am 2 and 9) ranges and limits. 



The analysis of solar neutrino 
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28| measurements involves considerable input from solar 



physics and the nuclear physics involved in the extensive chain of reactions that together are 
termed "Standard Solar Model" (SSM) [29j. Since the predicted flux of solar neutrinos from 



SSM is very well-established 



291 ] we know that the low energy p — p neutrinos are the most 



abundant and, since they arise from reactions that are responsible for most of the energy 
output of the sun, the predicted flux of these neutrinos is constrained very precisely (±2%) 
by the solar luminosity. The same is not true for higher energy neutrinos for which the 
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flux is less certain due to uncertainties in the nuclear and solar physics required to compute 
them. 

In fact, the true frontier for solar neutrino experiments is the real-time, spectral measure- 
ment of the flux of neutrinos below 0.4 MeV produced by p — p reactions. Measurements 
of the p — p flux to an accuracy comparable to the accuracy of the SSM calculation will 



significantly improve the precision of the mixing angle 20|, |2l[]. In addition, if the total 



p — p flux is well-know, measurement of the active component can help constrain a possi- 
ble sterile admixture [ll]. Therefore, an accurate phenomenological analysis for obtaining 
small modifications to the oscillation parameters, as we have illustrated in the Fig. HI can 
be really pertinent. In this context, the experimental challenge is to achieve low background 
at low energy threshold. In particular, a number of projects and proposals aiming to build 
p — p neutrino spectrometers are discussed in 32] where it has been shown that a large 



volume liquid organic scintillator detector with an energy resolution of 10 keV at 200 keV 
can be sensitive to solar p — p neutrinos, if operated at the target radiopurity levels for 
the Borexino detector, or the solar neutrino project of KamLAND 26| . In spite of the 
higher energy neutrinos being more accessible experimentally, the corrections to the wave 
packet formalism can be physically relevant for p — p neutrinos with energy distributed 
around an averaged value of E pa 10 — 100 keV. Following some standard procedure 36] 
for calculating the neutrino flux wave packet width a for p — p solar reactions, we obtain 
a = 10" 10 — 10~ 8 m = 0.5 — 50 (keV)' 1 . Such an interval sets a very particular range for the 
a parameter comprised by the interval 5-10" 5 — 0.2 for E pa 0.01 — 0.4 Me V which introduces 
the possibility for WP second-order corrections establish some not ignoble modifications for 
the p — p neutrino oscillation parameter limits. 

Turning back to the confront between the PW and the WP formalism, once we had 
precise values for the input parameters (P), L and E,we could determine the effectiveness 
of the first /second-order corrections in determining Am 2 for any class of neutrino oscillation 
experiment. For instance, the flux of atmospheric neutrinos produced by collisions of cosmic 
rays (which are mostly protons) with the upper atmosphere is measured by experiments 
)repared for obser ving v e <-> and v e <-* conversions. In particular, SuperKamiokande 



30] and MACRO 3JJ measurements also work on the <-> v T conversion. The neutrino 



energies range about from 0.1 GeV to 100 GeV which constrains the relevance of WP effects 



to an wave packet width a ~ 10 12 m. The accelerators experiments 33|, [34], [35] cover a 
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higher variety of neutrino flavor conversions where the neutrino energy flux stays around 
1 — 10 GeV the limitations to the wave packet width a (see the Table [I]) are analogous which 
makes the WP second-order corrections, as a first analysis, completely irrelevant. 

The most prominent contribution from the above discussion in determining the oscillation 
parameters and limits can come with the analysis of supernova neutrinos which, however, 
are not yet solidly established by the experimental data. Neutrinos from SN1987A in the 
Large Magellanic Cloud were detected by Kamiokande and 1MB detectors - only 19 events. 
Nowadays, SuperKamiokande, SNO, LVD, ICARUS, IceCube are expected to detect events 
from the next galactic supernova and improve the statistics, providing new information on 
neutrino properties and supernova. The main problem in studying neutrinos oscillation from 
supernova is the spectral and temporal evolution of the neutrino burst. In a supernova, the 
size of the wave packet is determined by the region of production (plasma), due to a process 
known as pressure broadening, which depends on the temperature, the plasma density, 
the mean free path of the neutrino producing particle and its mean termal velocity 36]. 
Neutrinos from supernova core with 100 MeV energy have a wave packet size varying from 
~ 5 • 10~ 16 m to ~ 10~ 14 m which leads to a wave packet parameter aE comprised by the 
interval 0.25 — 5 for which the second-order corrections can be relevant if aE > 1. There 
also are 10 MeV neutrinos, from neutrinosphere, with a wave packet size of approximately 
10" 12 m with aE « 50. 

In table [T] we summarize the results for the five well-known neutrino sources: reactor, 
solar, atmospheric, accelerator and supernova. We compare the size of the wave packet a 
determined by the region within which the production process of a neutrino is effectively 
localized by the physical process itself (for instance, from a cross section analysis) with 
the limits below which the WP second-order corrections set a more accurate value to the 
mass-squared difference Am 2 . It allow us to establish where/when the physical effects due 
to second-order corrections can be eventually be expected. 

The limits presented in table U represents, for instance, a diminution of 10%, 1% and 
0.1% from the value computed with the PW approach. However, a pertinent objection to 
the the representativeness of such values can be stated: it is important to observe that 
neutrino energy measurements cannot be performed very precisely so that this produces an 
effect competing with that of the finite size of the wave packet. If we set the energy un- 
certainty represented by 5E, the Heisenberg uncertainty relation states that 5E a ~ 1 and, 
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consequently, our approximation hypothesis leads to ~ ~ -~ -C 1. Realistically speak- 
ing, a typical neutrino-oscillation experiment searches for flavor conversions by means of an 
apparatus which, apart from the details inherent to the physical process, provides an indi- 
rect measurement of the neutrino energy (in each event) accompanied by an experimental 
error AE exv due to the "detector resolution". In case of AE j; xp < M ^ J_ the effective 

KJ -P E E aE 

role of the second-order corrections illustrated in this analysis can be relevant. On the 
contrary, AE exp > 5E demands for an averaged energy integration where the decoherence 
effect through imperfect neutrino energy measurements by far dominates. In a quantitative 
analytical analysis, this problem could be overcome by performing an analytical energy in- 
tegration which, in general, is not possible. In the present context, the current experimental 
values/measurements set some limitations to the applicability of our analysis by restricting 
it to the 7 Be and pe~p lines for solar neutrinos for which the energy flux is precisely defined 
^ AEex P ^ and, eventually, to some (next generation) reactor experiment and, certainly, 
to supernova neutrinos. 

As an additional remark, without comprehending the exact theoretical coupled to experi- 
mental procedure for determining the wave packet widths for a certain type of neutrino flux, 
we cannot arbitrarily assume, apart from the obvious criticisms to the PW approach, the 
modifications introduced by the WP treatment (particularly with second-order corrections) 
are irrelevant in the confront with a generic class of neutrino experimental data. Finally, 
from the phenomenological point of view, the general arguments presented in [11] continue to 
be valid, i. e. the above discussion has so far been limited to vacuum oscillations, where the 
mixing probabilities are described in terms of the mixing angle. In the solar neutrino case 
it is likely that interactions between the neutrinos and solar electrons affect the oscillation 
probability |37| . 



V. CONCLUSION 



We have reported about the intermediate WP prescription in the context of neutrino 
flavor oscillations from which, under the particular assumption of a sharply peaked mo- 
mentum distribution which sets an analytical approximation of order 0((aE)~ 3 ), we have 
re-obtained an explicit expression for the flavor conversion formula for (U)R and NR prop- 
agation regimes. By concentrating our arguments on quantifying the second-order effects of 
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such an approximation, we have observed that the existence of an additional time- dependent 
phase in the oscillating term of the flavor conversion formula and the modified spreading ef- 
fect can represent some minor but accurate modifications to the (ultra)relativistic oscillation 
probability formula which leads to important corrections to the phenomenological analysis 
for obtaining the neutrino oscillation parameter and limits. 

In particular, we have quantified such effects for determining the corrections to the 



Am 2 x sin 2 (26) curve for two reactor experiments 25|, |26j . The oscillation parameter range 
deviation from the PW values depends effectively on the product between the wave packet 
width a and the averaged energy flux E which characterizes the detection process. The im- 
portance of the second-order corrections which come from the WP construction can also be 
relevant in the framework of three-neutrino mixing. It is well diffused that the next question 
which can be approached experimentally is that of e3 mixing. A consequence of a non-zero 
U e3 matrix element will be a small appearance of u e in a bean of v^. for the particular case 
where Am 2 2 <C Am 2 3 (experimental data), and for E v ~ LAm 2 3 , ignoring matter effects, 
we can find |3i 



Pi^fj. —>■ v e ; L) = sin 2 (20is) sin 2 (6*23) sin 2 



Am 2 23 L 



(28) 



4E 

This expression illustrates that 6 13 manifests itself in the amplitude of an oscillation with 2-3 
like parameters. By assuming an intermediate wave packet analysis, fine-tuning corrections 
can eventually be relevant. 

Experimentally, since the modulation may be parts per thousand or smaller, one needs 
both good statistics and low background data. For instance the KamLAND experiment will 
significantly reduce the allowed region for Am 2 2 and sin (28 12 ) relative to the present results, 
where the second-order wave packet corrections can appear as an additional ingredient for 
accurately applying the phenomenological analysis. At the same time, the next major goal 
for the reactor neutrino program will be to attempt a measurement of sin 2 (26 13 ). It can be 
shown that the reactor experiments have the potential to determine 6 13 without ambiguity 
from CP violation or matter effects (by assuming the necessary statistical precision which 
requires large reactor power and large detector size). With reasonable systematic errors 
(< 1%) the sensitivity is supposed to reach about sin 2 (29 13 ) ~ 0.01 — 0.02 [3j and an 
accurate method of analysis, maybe in the wave packet framework, can be required. 

Turning back to the foundations for applying the intermediate wave packet formalism in 
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the neutrino oscillation problem, we know the necessity of a more sophisticated approach is 
required. It can involve a field-theoretical treatment. Derivations of the oscillation formula 
resorting to field-theoretical methods are not very popular. They are thought to be very 
complicated and the existing quantum field computations of the oscillation formula do not 
agree in all respects Q|. The Blasone and Vitiello (BV) model 0, Q to neutrino/particle 
mixing and oscillations seems to be the most distinguished trying to this aim. They have 
attempt to define a Fock space of weak eigenstates to derive a nonperturbative oscillation 



formula. Also with Dirac wave packets, the flavor conversion formula can be reproduced 381 

nn 

with the same mathematical structure as those obtained in the BV model [8J, H0[. In fact, 
both frameworks deserve a more careful investigation since the numerous conceptual diffi- 
culties hidden in the quantum oscillation phenomena still represent an intriguing challenge 
for physicists. 
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FIG. 1: The time-behavior of Osc(t) compared with the standard plane-wave oscillation given 
by cos for different propagation regimes. The additional phase @(t) changes the oscillating 

2 

character after some time of propagation. The minimal deviation occurs for j& ~ § which is 
represented by a solid line superposing the plane- wave case. We have used a E = 10 for this plot 
which was taken from reference [7|. 
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FIG. 2: We plot the behavior of the corrected phase A<&(£) = A£7t + 0(t) for different propagation 
regimes and we observe that the values assumed by 0(i) are effective only when the interference 
boundary function Bnd(£) does not vanish. By diminishing the value of the wave packet parameter 
aE (we also have used aE = 10 for this plot) the amortizing behavior is attenuated and the range 
of modifications introduced by the additional phase @(t) increases. This plot was taken from 
reference [7J. 
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FIG. 3: The time-dependence of the flavor conversion probability obtained with the introduction 
of second-order corrections in the series expansion of the energy for a strictly peaked momentum 
distribution (O(of)). By comparing with the PW predictions, depending on the propagation 
regime, the additional time-dependent phase A <!>(£) = AEt + G(t) produces a delay/advance 
in the local maxima of flavor detection. Phenomenologically, we shall demonstrate that such 
modifications allow us to quantify small corrections to the averaged values of neutrino oscillation 
parameters, i. e. the mixing-angle and the mass-squared difference. Essentially, it depends on the 
product of the wave packet width a by the averaged energy E. Here again we have used aE = 10. 
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FIG. 4: Neutrino oscillation parameter ranges excluded by two experiments: the Palo Verde reactor 
disappearance experiment (BOEHM 01) and the positive signal from the KamLAND collaboration 
(EGUCHI 03). We process the input parameters by following the WP treatment with first and 
second-order corrections (WP1 and WP2) where we have assumed the wave packet parameter a ~ 
0.342 which leads to a 10% correction to the value of Am 2 obtained with the PW approximation. 
In order to compare the results obtained with the WP approach with the PW prediction for the 
maximum excursion of the curve to the left, which corresponds to an important phenomenological 
boundary region, we have plotted the PW curve (dotted line) only for the smaller values of Am 2 
(since it has to be averaged on the energy for larger values). The parameter ranges simultaneously 
allowed by both experiments are represented by the not filled area. The second-order corrections 
introduce accurate modifications to the oscillation parameter ranges. The nearer is ex to 1, the 
more relevant is the contribution due to higher order terms in the Eq. (jlOp for determining such 
accurate limits. 
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FIG. 5: Relative corrections to the Am 2 value (obtained with the PW analysis) by assuming the 

WP treatments with first (WP1) and second-order (WP2) corrections in the energy expansion. We 

/ /p\ 

normalize the result by dividing Am 2 by bo g in ( 2 6>) where obviously (P) PW = (-P)wpi = (-P)wp2 = (P) 
is obtained from the experimental results. For instance, by following the WP treatment with 
second-order corrections, the correction to the phenomenological parameter Am 2 corresponds to 
a diminution of approximately 10% of the Am 2 PW value when a ps 0.342 (aE ps 3), 1% when 
a « 0.100 (aE ^ 10), and 0.1% when a w 0.032 (aE ^ 30). 
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TABLE I: Wave packet width a estimative for which the WP second-order corrections start to 
be relevant in obtaining an accurate value of Am 2 . We have computed the a limit by assuming 
10%(<r ~ 0.342), l%(a ~ 0.100) and 0.l%(a ~ 0.032), corrections to the value computed with the 
PW approximation. 
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